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1. Introduction 

The main goal of this paper is to solve the unsolved problem of the summation over Pomeron loops, in 
QCD. In proton proton collisions the main t-channel exchange is the BFKL Pomeron, which is a double 
t channel gluon exchange with "ladder" gluons in between, as shown in Fig. ||. The BFKL Pomeron can 
split into two daughter Pomerons and re-merge through the triple Pomeron vertex shown in Fig. |^, forming 
Pomeron loops. Due to the large size of the triple Pomeron vertex (see for example refs. ^) Pomeron 
loop diagrams give a significant contribution to the high energy scattering amplitude in proton proton 
collisions, which is comparable to the amplitude of the basic diagram of Fig. |l]. As such, the scattering am- 
plitude of hadronic reactions requires an accurate estimate for the summation over Pomeron loop diagrams 
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Figure 1: The BFKL Pomeron structure. 



Figure 2: The triple Pomeron vertex. 



to be taken into account (for an example of such diagrams, see Fig. ^, Fig. 14 and Fig. 15). Hence the 
correct algorithm for the summation over Pomeron loops is a result which is in high demand, not just from 
a theoretical perspective, but also from an experimental point of view. The main practical application of 
this result is in exclusive diffractive Higgs production, at the LHC. 
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Figure 3: Diffractive production of the Higgs boson through double t-channel gluon exchange. On the left is exclusive 
Higgs production with large rapidity gaps (LRG) between the Higgs and the emerging protons. On the right is the 
production of Higgs with extra production which spoils the LRGs, arising from additional inelastic scattering. 

The keenly awaited result of the discovery of the Higgs boson at the LHC, is expected to emerge from 
smashing together two protons, in so-called diffractive Higgs production. Unfortunately, proton-proton 
scattering results in the production of many other unwanted particles. This makes the process of detecting 
the Higgs very problematic. The desired result is the production of the Higgs with two large rapidity gaps 
(LRG) between the Higgs and the emerging protons as shown in the left hand diagram of Fig. 0, which 
ensures that there is no additional production. The survival of these two large rapidity gaps is quantified 
by the survival probability. The difficulty in isolating the Higgs signal is characterized by the small value 
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of the survival probability, which has so far been estimated to be small, and could be even less than 1 % 
(see detailed estimates in refs. fs], |^). 

Due to the size of the strong coupling, the dominant mechanism for diffractive Higgs production is 
through the exchange of a t channel gluon between the scattering protons. To ensure that there is no 
further production, a second t - channel gluon exchanged between the scattering protons is needed to 
cancel the color flow, as shown in Fig. |3|. This colorless double gluon exchange cancels the possibility of 
the production of additional particles. 

The double gluon exchange evolves to the "ladder" structure, as gluons are exchanged between the 
two t channel gluons, forming the rungs of the ladder (see Fig. |^. This structure is the so called BFKL 
Pomeron. The energy levels of the BFKL Pomeron are labeled by the BFKL eigenfunction to (n, z^), where 
n is an integer and i/ is a conformal variable which one integrates over, when calculating Feynman diagrams. 

Unfortunately, extra parton showers between the scattering protons is inevitable, leading to the pro- 
duction of additional particles that fill up the large rapidity gaps (see Fig. ^ right). Hard re-scattering 
corrections, namely Pomeron branching which forms loop diagrams and fan diagrams of the type shown 
in Fig. 1^ also contribute to the problem of additional unwanted production. Pomeron loop diagrams and 
fan diagrams of this kind are called "enhanced" diagrams. The suppression factor needed, to remove the 
effect of enhanced diagrams on the Higgs cross section, is the "enhanced survival probability" < {S"^^^] >. 

Pomeron loop diagrams are extremely useful because they include the contribution of "fan diagrams" , 
Pomeron loops and the contribution of non interacting Pomerons as shown in Fig. ^. As such, summing 
over the complete set of Pomeron loop diagrams will provide the suppression factor needed for screening 
out all unwanted particle production, which stems from Pomeron enhanced diagrams. In other words, the 
summation of Pomeron loops provides the full set of hard re-scattering corrections to the enhanced survival 
probability < \Sl^^\ >. 

In ref. Q the one loop Pomeron diagram was calculated in QCD. The result was used to fix the pa- 
rameters of the mean field approximation (MFA) of Mueller, Patel, Salam and lancu (MPSI) (see ref. [^), 
and using this the summation over Pomeron loops was estimated using the toy model of Kovchegov Q . In 
this way the survival probability was found to be potentially as low as 0.4%. This implied that Pomeron 
loop diagrams give a significant contribution to the survival probability. In ref. |0] we derived a useful 
expression for the triple Pomeron vertex in the momentum representation, which is the crucial ingredient 
for calculating any diagram with Pomeron branching and Pomeron loops. Using this we showed that a 
diagram with an arbitrary number of Pomeron loops reduces to the diagram of non interacting Pomerons. 
This important step forward, outlined the first method for the summation over Pomeron loop diagrams 
completely in the QCD approach, and not in the mean field approximation. 
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Figure 4: The first Pomeron loop diagram gives the simultaneous contribution of the fan diagram, the loop diagram 
and the contribution from non interacting Pomerons. 

In a recent publication by M. Braun, the Pomeron loop diagram was calculated using the original 
exact form of the triple Pomeron vertex, which was first derived by Korchemsky Q and at the same time 
by Bialas, Navelet and Peschanski (BNP) in ref. It was suggested in ref. that use of the Korchemsky 
expression for the triple Pomeron vertex, leads to a more accurate result for the Pomeron loop diagram. 
Therefore the approach suggested in ref. |^ has been adopted in this paper. 

The approach used to find the Pomeron loop amplitude, can be extended to the calculation of more 
complicated diagrams (see for example Fig. ^l]), with an arbitrary number of Pomeron loops in the QCD 
approach instead of the MFA formalism. Using the new algorithm developed in this paper, the amplitude 
for multiple loop diagrams derives from an iterative expansion of the simple Pomeron loop amplitude. In 
this framework, a general formula can be found for the amplitude of the diagram with an arbitrary number 
of Pomeron loops in QCD, which is a function of the number of loops in the diagram, and the rapidity 
gap between the scattering protons. If the production of the Higgs boson is included in the diagram, the 
formula is also a function of the mass of the Higgs boson. This provides the mechanism for summing over 
Pomeron loop diagrams in the precise QCD formalism, instead of the toy model approach of ref. ||5|, ^. In 
this approach it was found that the main contribution to the amplitude of the general enhanced diagram, 
is equivalent to the amplitude of the diagram with non interacting Pomerons, with renormalized Pomeron 
vertices. 

This is a considerable step forward theoretically, since it makes possible the summation over Pomeron 
loop diagrams in QCD to any order. Experimentally this means that an accurate prediction for the 
exclusive cross section for diffractive Higgs production is possible that includes the suppression factor 
needed to screen out the full set of re-scattering corrections, namely the enhanced survival probability. In 
this new formalism it was found that the contribution of enhanced diagrams to the survival probability is 
substantial, confirming the results found in refs. |^, The results show also that the survival probability 
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is very sensitive to the choice of the strong couphng, and as such it decreases as ct^ increases, in agreement 
withref. |]. 

This paper is organized in the following way. In section ^ the notation and conventions used through- 
out this paper are listed. The amplitude for the basic diagram of Fig. 0, which shows diffractive Higgs 
production through single t-channel Pomeron exchange is derived. Section |^ is devoted to the derivation 
of the amplitude of the Pomeron loop, shown in Fig. |^. A pedagogical approach is taken to explain the 
style of integration over the conformal variables in the loop. This forms the basis of the techniques used for 
more complicated diagrams in later sections. In section Q the main idea used to derive the formula for the 
general multiple Pomeron loop diagram is outlined, which is an iterative technique. Using this iterative 
approach, in section |5|, more complicated multiple Pomeron loop diagrams, such as Fig. 14, Fig. ^ and 
Fig. ^ are calculated. Finally in section ^ the general expression for the amplitude of the diagram with an 
arbitrary number of Pomeron loops is derived. This provides the tools necessary to sum over the complete 
set of Pomeron loop diagrams up to any order. The summation over enhanced diagrams forms the basis 
for section |^, to estimate the enhanced survival probability < | S'^^y^ \ > . These results are discussed in the 
conclusion in section |8[ In the appendix section a brief overview of the derivation of the triple Pomeron 
vertex is given, mostly following the strategy of Korchemsky in ref. ffl, and the amplitude of the first fan 



diagram of Fig. 26 is calculated. 



2. The Pomeron propagator 

The following conventions and notations will be used, some of which are based on the paper of ref. 1^ by 
M. Braun. The Pomeron propagator shown in Fig. ^, which is a t channel exchange between two color 
dipoles in QCD, with rapidity values y and y' , is denoted by the expression; 
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Figure 5: Coordinates for the two reggeized gluons in the BFKL Pomeron propagator, in color dipole scattering. 
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where as shown in Fig. ||, ri and r2 denote the initial and final coordinates of the two reggeized gluons, 
and R is the center of mass coordinate of the two reggeized gluons. The conformally invariant expression, 
for the coupling of the BFKL Pomeron to the QCD color dipole is; 



where rij = — rj and rjo = ri — R. The confer mal weights are given by the expression 



1 + n 

7 = —^+iv] 7 = 1-7*; n G EM (2.3) 

Two important properties of the conformal weights are 



7— 7 = n 7 + 7 = 1 + 2iv (2-4) 

where n is the conformal spin, and v is the scaling dimension of the state. Wherever 7 appears, it 
is intended to be a shorthand for the set of two numbers {n , z^}. The above expression for the Pomeron 



propagator in Eq. (2T) can also be written in the representation of complex angular momentum j = l+to 



instead of rapidity y by the following Mellin transformation; 

where in u; representation, the Pomeron propagator is given by the following expression; 

= jd^RjcfR' jv^E^EL^Qu^ij) (2-6) 
where P7 is a shorthand notation for 

j V-f = j djhi-f); where /i(7) = A 7 - i ^ (2.7) 

The contour shown in Fig. ^ consists of the imaginary 7 axis from ±ioo, and the semi circle at 
infinity, to the left of the imaginary 7 axis. Cy encloses all singularities in the integrand of Eq. (p.6|). The 
integrand vanishes on the semi circle at infinity, such that it is sufficient just to replace 
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It is more economical to calculate diagrams with BFKL 
Pomeron states in terms of v, instead 7. For the integration 
limits e + ioo < 7 < e + i 00 (as e ^ 0), Eq. ( |2.3D gives the 
corresponding limits of integration for the variable i^, as 
— 00 < < 00, and one should sum over all real positive 
integers n. In this way in Eq. (2.6), one can replace the 



integration over 7 with the integration over v using the 
notation 



Figure 6: Contour enclosing singularities for the 
integration over the conformal variable 7. 



d'y 



E 



du 



(2.9) 



The notation Yl'n=~oo I^oo corresponds to the integration over the quantum numbers associated 
with the continuous unitary variable irreducible representations of SL ( 2 , C ) , defined in Eq. ( ^.3[) (see 
refs. m, |ll|] for a detailed explanation). Thus in terms of u; 
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dv h{n , I') 



(2.10) 



where h{n , u) 

The conformal propagator (7^ (7) takes the form; 
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(2.11) 
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where 



\{n , v) = — X 



uj — uj [n ,v) 
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16 h{n + 1 , u) h{n — 1 , v) 
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16 + + 1)2 /4| + _ 1)2 



(2.12) 



(2.13) 



a; (n , z^) are the eigenfunctions of the BFKL equation, which represent the energy levels of the BFKL 
Pomeron, and these are given by the expression 
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Figure 7: Diffractive production of the Higgs boson through t-channel Pomeron exchange. The process Pomeron + 
Pomeron — > Higgs proceeds mostly through the top quark triangle. 



Loin,iy) = a,{V'(l)-3f?eV(7)} = «a2V'(l)-V'(7)-V'(l-7)} (2.14) 
where the function 



is the Di-gamma function, and where T (x) = (x — 1)! is the Euler Gamma function. Throughout 
this paper, the notation 



a. 



(2.16) 



where Nc is the number of colors in the SU {Nc) color group, or abstractly speaking, Nc is the number 
of generators for the standard representation of SU {Nc). For QCD, Nc = 3. Here, Og is the strong 
coupling for QCD interactions, and it depends inversely on energy. Wherever a value for as is used in a 
calculation, its value for a specific choice of energy scale will be stated explicitly. 



Substituting for the conformal propagator Eq. ( 2.12| ) in Eq. (p.5|), and using the integration measure of 
Eq. ( ^.10| ) to integrate over the conformal variable, the amplitude for diffractive Higgs production through 
the t channel Pomeron exchange process shown in Fig. |^ is; 
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A^o){^,5yH\Fig.^) =^Y1 / di^ h{n,u)X{n,u)eicp{u;{n,u){A-6yH)}E^E'_yAH (2.17) 

where A = y — y' is the rapidity gap between the scattering dipoles shown in Fig. 0, and Syn = 
In (M^/4so) is the window of rapidity occupied by the heavy Higgs boson {Mh is the mass of the Higgs 
boson and so = IGeV^). The notation ^(q) represents the fact that there are = generations of 
Pomeron branching, since later on the notation ^(a?) is used to label the diagram with N generations of 
Pomeron branching. The numerical factor of in front of Eq. ( |2.17 ) takes into account the couplings 



of the two reggeized gluons in the BFKL Pomeron structure in Fig. to the color dipoles. Since there are 
two gluons which each couple twice at both ends to the QCD color dipole, this leads to a factor which is 
proportional to a^. The factor of 1/4 accounts for the 4 degenerate diagrams in Fig. ^, which arise due to 
the different ways of coupling the reggeized gluons to the color dipoles. Hence one divides by a symmetry 
factor which ensures that identical diagrams are only counted once. A more detailed explanation of this 
coefficient can be found in refs. ^ |l^, 12]. Ah in Eq. ( 2.17| ) represents the amplitude for the subprocess 



where Pomeron + Pomeron Higgs, which is given by the expression 14, 15, 16, O, Q 



3 TT 

= 6.89 X 10-^ GeV-2 (2.18) 

where Gi;' is the Fermi constant. Here a typical value for Og, at the scale of Mz, the mass of the Z 
boson, is used. It is expected that the Higgs will be produced with a mass of approximately lOOGeV/c^, 
which would give a value for the strong coupling constant ~ 0.12. This corresponds to a Z particle mass 
[|19|, of Mz = 90.8 ± 0.6 GeV/c^. The leading order contribut ion to this process is where one of the two 
reggeized gluons in the BFKL Pomeron states produces the Higgs boson through the quark triangle shown 
in Fig. ^, for which the top flavour dominates. This is on account of the quark + quark — > Higgs vertex 
in Fig. ^ which is proportional to the mass of the quark flavour, so the top flavour which is the heaviest 
("^top = 175GeV/c^) contributes the most. 

The BFKL eigenfunction specified in Eq. ( p. 14 ) decreases sharply as n decreases. In fact, the BFKL 



eigenfunction remains positive at high energy, only when n = 0. For this reason, solutions with n ^ are 
negligible, and from here onwards will be ignored. Note that from the definition of Eq. ( p. 14 ) in the case 
when n = 0; 

uj{u) = as ^^2^l;{l)-^l^(^^ + iu^ -ijl^^-iu^"^ (2.19) 

which has a saddle point at u = 0. Hence for the BFKL eigenfunction one can use the expansion 
around the point = 
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(2.20) 



so that Eq. (2.17) simplifies to 



^(o)(A,^,.|Fig.0) =^,j_J^-^^. 



X exp I uj (0) (A — dyu) — -^^"^ ^" {^) (^ — ^Uh) + \n[EE') \ Ah (2.21) 
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(2.22) 



The integration over v in Eq. (2.21) can be solved by using the method of steepest descents. In this tech- 
nique the exponential is expanded around the saddle point Vsp (in this case Vsp = ln{E E') /u" (0) (A — dyn)) 
and the remaining part of the integrand is fixed by setting the integration variable v = Usp- In this way 
the integration reduces to a Gaussian type integration over v, and using the result that dx e~"^^ = 
7ri/V2 X a'^''^ (« / 0), Eq. (|2^ ) yields the result; 



A(0) (A,dyH|Fig. I) 



47r4 
X ( 1 + 



Lo"{Q){^-5yH] 

\n{EE') 
c^"(0) {^-Syu) 



3/2 



exp <^ u (0) (A - Syn) 



ln2 {EE'] 



2cu"(0) (A-5yH) 



(2.23) 



In the case where the v saddle point Ugp = \n{E E') /{uj"(0) (A — Syn)) is small, then one can recast 
Eq. ( |2.23| ) in the following asymptotic form, namely 



^(0) (A,5y/^|Fig. 0) 



(2vr) 



H 



27r2Ar2 



(u;"(0) {^-5yH)f'^ 



1.57 X 10"® GeV"2 (a^ = 0.12) 
2.17 X 10"' GeV-2 («, = 0.2) 



(2.24) 



The following values were used; 



a;(0) = AoLs In 2; u;"(0) 
A =19; 5yH 



28a,C(3); 
In 



Ml 



4so 
7.824 



(2.25) 
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Figure 8: Plot of A(o) ( A, (5^^ |Fig. |7|) derived in Eq. ( ^.24D 



against the rapidity gap A between the scattering dipoles. The 
values for the rapidity gap A go from Syn = In (M^/4so) for 



the Higgs mass Mh = 100 GeV and sq 
typical LHC rapidity A = 19. 



IGeV , up to the 



where C (n) is the Riemann-zeta function, 
and where the rapidity gap A between the scat- 
tering dipoles is based on the energy y/s = 14 
TeV, which is typical for proton collisions at 



the LHC. The values given in Eq. (2.25) will be 
assumed throughout this paper. The difference 



in the values of Eq. ( 2.24 ) shows that the am- 
plitude depends critically on the choice of the 
strong coupling. The plot of the single Pomeron 
amplitude derived in Eq. (|2.24|) against rapidity 
is shown in Fig. for Ug = 0.2. The plot shows 
that the amplitude of Fig. |^ is very sensitive to 
increases in the rapidity gap A. 



3. The simple Pomeron loop 

In this section, the Feynman amplitude for the diffractive production of the Higgs boson, from the Pomeron 
loop diagram is derived. The diagram shown in Fig. is the scattering of a color dipole at rapidity y off a 
color dipole at rapidity y ', with the exchange of a BFKL Pomeron state, where the correction of one loop 
is included. Its amplitude is denoted (u^lo' ,uji), since there is one generation of Pomeron splitting, 
and likewise one generation of recombining of two Pomerons into one, to form the loop. 



^(1) (u;,u;',c^i|Fig.|) =^jv^j V^' E'^^,g^ {^) g^, {^') D(^^ {^,^'\u:,uJi) Ah (3.1) 

The contribution Df^i^ (7, 7'|ti;,ti;i) of the Pomeron loop to the amplitude of Fig. |9|, is written in uj 
representation as a function of the two conformal weights 7 and 7 ' of the Pomeron propagators attached 
to the loop, namely; 

(7,7^,(^1) = J cfRj (fR' d^,){-f,Y\u;,uJi) (3.2) 

where R and R ' are the center of mass coordinates of the two Pomeron propagators on either side of 
the loop, and 
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Figure 9: The Pomeron loop diagram. 



d{i) / (l^^ij d'^R2j T>-fijv-f2gu;i (71) 9io-ioi (72) 

xV iR,Ri,R2h ^ ^i,^2)V {R',Ri,R2hi,^2 ^ 7') (3.3) 

5(1) in the denominator is a symmetry factor which prevents counting identical diagrams more than 
once, and its value is determined below. V {R, -R2I7 7i) 72) is the vertex for the splitting of the BFKL 
Pomeron state labeled 7}, into the two BFKL Pomeron states labeled by {i?i,7i} and {-^2)72}) and 
it is known as the triple Pomeron vertex, shown in Fig. ( 7 = conformal variable and R = center of 
mass coordinate). 

The total order of the symmetry group of Fig. |^ is S'(i) = 16. This can be seen by observing that for 
both vertices of the loop, there are 2 permutations of the Pomeron lines which form the 2 branches of the 
loop. In addition for each pair of Pomeron lines there are 2 permutations formed by swapping the upper 
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with the lower vertex. This makes a total symmetry factor of 4. The remaining 4-fold degeneracy stems 
from the triple Pomeron vertex shown in Fig. |l^, whereby swapping the reggeized gluon lines (5) and (6) 
leads to an identical diagram. However swapping lines (3) and (4) leads to two different diagrams called 
the planar and non -planar diagrams shown in Fig. |l^ and Fig. [l^. These two diagrams are non identical, 
and give different contributions, where the latter is suppressed by 




Korchemsky pointed out in ref. 
that by the condition of conformal invari- 
ance, the triple Pomeron vertex can be 
factorized into two pieces which depend 
only on the center of mass coordinates, 
and the conformal variables respectively, 
in the following way; 



Figure 10: The triple Pomeron vertex. 



V{R,RuR2h^ 11,12) = Roi''''Ru''''R2o''''Roi^''Ru^''R*20^''^il\li,l2) (3.4) 

where R^i = R — Ri {i = 1,2) and R12 = Ri — R2, and where for example Aoi = 7 + 71—72! and 
A12 = 71+72 — 7) with an equivalent definition for the A in terms of the 7. It was found in refs. [||, 
that; 



r(7l7i,72) = 167(1-7)7(1-7) 



TT 



X <{ S7 (7 I 71 , 72 ) + A (7 I 71 , 72 ) (x (71) + X (72) - X (7)) )■ (3.5) 



planar piece Fig. |l^ non-planar piece Fig. ^ 



where x (7) = {ip (1) - tp {j)} (3.6) 

where ijj (7) is the di-gamma function. The triple Pomeron vertex comes in two types, namely the 
planar and the non planar diagram shown in Fig. |l^ and Fig. |l^ respectively. The first term in curly 
brackets is the contribution of the planar diagram, and the second term is the contribution of the non 
planar diagram. The explicit expressions for $7 (7 | 71 , 72 ) and A (7 | 71 , 72 ) and their derivation, can be 
found in the appendix. 
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Figure 11: Planar diagram [Q. 



Figure 12: Non planar diagram 



By inserting Eq. { \iA\ ) into Eq. (3.3) and evaluating all the integrals over the center of mass coordinates, 



it was shown in ref. M that this leads to the following condition; 



d(i) {7 , ^'\uJ,uJi) = d(i){j\iv,ivi) j-^Sn,n'S {i^ - 1^') S^^^ {R- R') (3.7) 
where (7 |w, tJi) = ^ /2?7i (71) 5<.^-a;i (72) T (7I71, 72) T (7I71, 72) (3.8) 



Now inserting the loop amplitude of Eq. ( |3.8D into the scattering amplitude of Eq. ( |3.1| ) , the integration 
over 7' can be eliminated by virtue of the Dirac delta function 5(7 — 7')- Hence after introducing the 



explicit expressions for the integration measure and the conformal propagator given by Eq. (2.10) and 
Eq. (|2.12|) respectively; 
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r (7I71, 72) r (7171,72) (3.' 



The amplitude of the Pomeron loop diagram can be re-expressed as a function of the rapidity instead 
of the angular momentum j = \ + u), using the Mellin transform. The rapidity gap filled up by the loop 
is Ai = yi — y[, where yi and y[ are respectively the rapidity values of the upper and lower vertices of 
the loop (see Fig. ^ . In the case of inelastic scattering which results in the production of the heavy Higgs 
boson from one branch of the loop, there will be a window of rapidity 6yH which the Higgs boson occupies. 
Therefore dyn should be the minimum rapidity gap which the loop occupies, such that the energy of the 
scattering is enough to produce the Higgs boson. This affects the upper and lower limits of the rapidity 
variables yi and and hence the choice of integration limits are as they appear below in Eq. ( p^ ; 
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f^^^iy^y.) [^^{yi-y') [^My.-yi) A,,. (u:,u:[u;, 
J 2-Ki J 2m J 2m ^ ' ' 

fjv^X' (7) e'^W^d(i) {7\A,SyH) E,E'_^Ah 



Fie 



(3.10) 



where A = y — y' \s the rapidity gap between the scattering dipoles. The loop amplitude in rapidity 
representation is given by the expression; 



^^(1) {i\^,^yH) 



1 

16 

y 



D71 j P72 A (71) A (72) r (7 1 71 , 72) r (7 1 71 , 72) 

yi-&yH 



XI dyi I dyi exp{(a;(7i) +^(72) - ^(7)) (yi -2/1') -a;(7i)(5y//} (3.11) 
'y'+^yn Jy' 

Assuming that the leading contribution at high energy stems from the region where the conformal spins 
n = ni = 77-2 = 0, and introducing the exphcit expressions given in Eq. ( 2.10| ), Eq. ( 2.11D and Eq. ( |2.13D 
for the integration measure and the conformal propagator, the scattering amphtude of Eq. ( p. 101 ) reduces 
to 



^(i)(A,5yH|Fig.i =-^/ duh{u)X^{u)e'^^^^'^d^,){i.\A,SyH)E,EL,AH (3.12) 
where (when the conformal spins ni = 712 = 0) the loop amplitude simplifies to; 



dii)ii^\A,5yH) = ^ 



00 roo ^2 ^2 

dui / dl'2 7j 7j 

-00 J -00 {ul + 1/4) ^ (i/| + 1/4) ^ 



T{v\vi, 1/2) 



ry ryi-Syn 

XI dyi I dylexp{{uj {I'l) +UJ {1^2) - (i^)) {yi - yi) - i'^i)SyH} (3.13) 



y'+Syn Jy' 



A convention used in Eq. ( |3.13| ) and throughout this paper, is that wherever u appears alone without 

n, it has been assumed that n = and it has been suppressed. For example F (7I71, 72)„^„^^„2=o 

is labeled F (z^|z^i, 2^2)1 and w (n = 0, i^) is labeled w(i^). In the region where n = ni = 77-2 = 0, then 

r (7l7ii72)„=„^=„2=0 = r (-;/| -1^1,-1^2), and hence in Eq. ( |3.13| ), the notation 
2 

F {uiui, V2) = r (z^|i^i, 1^2) r {v\vi,i'2y = r (z^|z^i, 1^2) f (-v\ - ui, -1^2)- 
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Two main regions are considered when evaluating the {v,vi^i'2} integrals of Eq. ( 3.12| ) and Eq. ( ^.131 ). 
Region I is {z^, \vi\, |z/2|} = {0, 1/2, 1/2} and region II is z/i, U2} = {1/2, 0, 0}. To find the contribution 
from region I, after integrating over the rapidity variables yi and y/, the loop amplitude of Eq. ( 3.13 ) takes 
the form; 



210 



X e 



7r° 



.l^{l^2) SlJH 



00 1/ 



f + l/4)'{ui + l/Af 



T{iy\i^i, 1/2) 



X {A-6yH + P - Pe 



(3.14) 




Figure 13: The integration contour C (C) 
which is closed on the upper (lower) half plane, 
in order to enclose the pole ii/ — 1/2 (—1/2). 



where in Eq. (3.14) the identity l/x" = l/r(n) x 
e~^^ f3^~^ has been used. The singularities in the 
integrand at { | z^i | , 1 1^2 1 } = 1/2 suggest closing the ui and 
1^2 integration paths along the contour C (or C) shown in 
Fig. 13 which runs along the real axis from —00 to +00 and 
along the semi circle at infinity in the upper (lower) half 
plane. C (C) encloses the points {iz/i,zz/2} = 1/2 (—1/2), 
and the two solutions which stem from both contours are 
identical. Therefore the contour C is chosen, so that the so- 
lution is 2 X (27ri)^ (residues at {ivi,ii'2] = 1/2), where 
the extra factor of 2 takes into account the identical con- 
tribution from the residues at {11^1,11^2 = ~l/2}- The non- 
singular part of the integrand of Eq. (3.14), namely 
(1/2 + (1/2 + tends to zero as l/i^^f|, hence 



the convergence on the semi-circle at infinity is satisfied, al- 
lowing C to be used as the integration contour. This choice 
of contour will be used repeatedly in similar calculations 
throughout this paper, therefore wherever a specified integration contour is labeled C, it is assumed that 
the path shown in Fig. 13 is intended, and that the non-singular part of the integrand has a good conver- 
gence on C. 



The triple Pomeron vertex in region I which should be inserted into Eq. ( |3.14| ) is derived in the appendix 
in Eq. (A. 22). After taking into account all the singularities that stem from Eq. ( |A.22| ) at {ii'i,ii'2} = 1/2, 

and the singularities that stem from Lij (zzvi^2 — * 1/2) = ols/ (1/2 — ii/1^2) and x (^^^1,2 1/2) = 1/ (1/2 — 
(see Eq. ( ^.19| ) and Eq. (|3.6| )), the integration over vi in Eq. ( ^.14 ) takes the form; 
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(1/2 + ,.A/2 - ..0' {w^A 

where the derivative acts on the exponential to bring down a factor of (1/2 — ivi)~'^ . The (1/2 — iui)^"^ 
term remaining is canceled by the Jacobian of the coordinate transformation u = ckg (1/2 — ivi)^^ + 
q:<j (1/2 — ^2^2)""*^ —"^ (^'2) 5y_ff//3, for which the Jacobian is |0u/9z^i| = {1/2 — ivif' / ckg- Instead integrating 
over u yields an integration of the form §(jdu exp {—iu (3) = irid (/3) or f^du exp {—iu (/? — A + Syn)) = 
2TTi6 (/3 — A + dyn)- After taking the residue at ivi = 1/2, the integration over 1/2 has now reduced to; 



{l-ivi-W2-'2,/Nl) ty^-.i/2 r .( 1 2/iV, 



^ dV2 ^-^5 '- V- "~ i' dV2 f (1^2) 1 n T 

Ic (1/2 + iu2f (1/2 - 11^2)^ Jc V (1/2 - ii^2f (1/2 - iu2)^ 

= 2m t}lfi2) = 1/2) _ ^^/(^) {^U2 = 1/2)^ 

Overall using the above described method for taking the integrals over vi and 1^2, and evaluating the 
integration over /3, the loop amplitude of Eq. (|3.14 ) becomes; 



dl^{u\A,6yH) =d(^^ + zu^ i}''"") x(^){^(^^) + ^^'(^)(A-5yj^)|e-M'^^« (3.17) 
where d = (l - ^) (l - ^) (3.18) 



Now inserting the result for the loop amplitude of Eq. ( [3.17 ) into the scattering amplitude of Eq. ( 3.12| ), 



the steepest descent method described in section |2| can be used for the v integration in Eq. ( |3.12D , so 

that overall the contribution to the scattering amplitude which stems from the region {i/, |z^i|, I1/2I} = 
{0,1/2,1/2} is; 

(2ttV^'^ a"^ dAu ec^(o) (A-5s,h) fin 1 

= 2.34 X 10"^^ GeV"2 (q^ = 0.12) 

1.3 X 10"^ GeV-2 ^ 0.2) (3.20) 
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The results of Eq. ( |3.2[l| ) show that the Pomeron loop in region I is small compared with the basic 
diagram of Fig. ^ (see Eq. (|2.24| )). It turns out that the dominant contribution to the Pomeron loop 
amplitude stems from region II, namely {|z^|, z^i, 1^2} = {1/2,0,0}. In this region the expansion ofEq. ( g30D 



for the BFKL eigenfunctions uj (vi) and uj {U2) in Eq. ( |3.13| ) can be used. The triple Pomeron vertex of 
Eq. ( [A. 231 ) derived in the appendix for region II should be inserted into Eq. ( 3.13| ). The expression of 



Eq. ( |A.23| ) is novel for the following reason. From Eq. (^]^) the triple Pomeron vertex in region II takes 
the following form Q (n = ni = n2 = at high energy) ; 



r(|i/||zyi,i/2) (region II) = F ( -|0,0 ) = 16 ' "'^^^ 



2' \ TT 



2 



ti/- 



1 .\Vl . V ( ^ fl. \ 27r./l, \ f f . 1 



planar piece Fig. |ll] non-planar piece Fig. |l2| 



O(i|0,0) and A(i|0,0) were calculated in the appendix (see Eq. ( |A.11| ) and Eq. ( |A.15| )) and they 
both have a second order pole at iv = 1/2, canceled by the (1/2 — iv)^ term in front in Eq. ( |3.2lD . However 
an extra pole from the non planar piece arises due to the singularity which stems from x (^^ ~^ 1/2) 
(1/2 — (see Eq. (p.6|)). Therefore overall for region II the divergent part of the triple Pomeron vertex 

stems from the non planar piece, and the contribution from the planar piece is non-singular. Hence a 
remarkable feature of the Pomeron loop arises which has not been taken into account before, namely that 
the dominant contribution to the Pomeron loop amplitude comes from the non planar piece A (z^|z^i, 1/2) 
of the triple Pomeron vertex shown in Fig. |l^. This property has been so far ignored because previous 
publications assumed that ^ 00, and since in Eq. dl^) the non planar piece is suppressed by 2tt/N^ 



compared to the planar piece, it was neglected. Taking this into account and after evaluating the z^i 
and 1J2 integrations in Eq. ( 3.13| ) using the method of steepest descents, the contribution of region II 



({|i/|, vi, 1/2} = {1/2, 0, 0}) to the loop of Fig. ^ simplifies to; 



dll.{u\A,6yH) =j-j——rj-j:: r^/ dyi dy[—j- — (3.22) 

^ ^ 0-/2 + {1/2 - iv) Jy,+SyH Jy' (Ai - 5yi/)^/2 

2^0^ 



Inserting Eq. ( p.22| ) into Eq. ( 3.12| ), taking into account the singularities that stem from uo [iu — > 1/2) 



(1/2 — iv) ^ (see Eq. ( ^.19)) ) and closing the contour of integration on the path C shown in Fig. 13, the v 
integration takes the form; 
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, 1/2 [a,(A-Ai) 



c 



[1/2 + iuf {1/2- iuf I 1/2 -iu 



. dv = ^ — exp <^ ^ \ (3.24 

c {1/2 + wf {1/2- ivf \as d/\ J 1/2-tu j 

The remaining (1/2 — iv)~'^ term is canceled by the Jacobian of the transformation u = as/ (1/2 — iv), 
and integrating over u yields the derivative of the Dirac delta function 27ri5^ (A — Ai) /a^, which is ab- 



sorbed by the integration over the rapidity variables {yi,y[] in Eq. ( 3.12 ). Overall the contribution of 



region II ({|z/|, i^i, ^2} = {1/2,0,0}) to the scattering amplitude of Fig. is; 



^ a aAu f-1 { e^^^^HA-^J/H/s) 1 

,^,,,„|pi, i (A - .„) {-1±) ]^.^^^-—^^ ,3.25) 

= 2.05 X 10"^° GeV-2 (a, = 0.12) 

2.14 X 10""^ GeV-2 (a, = 0.2) (3.26) 

An obvious divergence in the integrand of Eq. (|3.1l| ) originates from the region where the conformal 
weights @ ni = it 1 and n2 = ± 1 (see Eq. (^l3|)). However it was proven in ref. Q that BFKL Pomeron 
states with odd ni or n2 cannot couple to BFKL states where n is even, and therefore this particular 
divergence does not arise. 

The complete scattering amplitude for Fig. ^ is given by the sum of the two contributions of Eq. ( p. 20] ) 
and Eq. ( 3.26| ), namely; 



(A, <5y^^|Fig. I) = ^(\) (A, 5yH\Vig. i + (A, |Fig. |) (3.27) 

= 2.28 X 10"^° GeV-2 (a, = 0.12) 

2.15 X 10"^ GeV"2 (a, = 0.2) (3.28) 

Comparing the results of Eq. ( ^.201) and Eq. ( |3.26| ), it is clear that for the simple Pomeron loop 
amplitude of Fig. |9|, the dominant contribution stems from region II. Prom an observation of Eq. ( ^.26 ), 
the simple loop amplitude can be written as 



All. (A, JynlFig. i = A;(i)e2-(0)A (3.29) 
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where contains ah the other terms included in the simple loop amplitude. From Eq. ( |3.29| ) it is 
clear that the dominant contribution to the loop amplitude is equivalent to the amplitude of 2 non inter- 
acting Pomerons. This can be seen from Fig. |9|, where by taking the 2 branches of the loop outside, one 
observes the exchange of 2 non interacting Pomerons. This was first noted by A.Mueller, who commented 
that at high rapidities, loop diagrams reduce to the exchange of multiple non interacting Pomerons, with 
renormalized Pomeron vertices (see ref. pi). 



Comparing the results of Eq. (|2.24 ) and Eq. ( |3.28D , the Pomeron loop of Fig. |9| is potentially very 
close to the basic diagram of Fig. |^ in amplitude, and suggests that Pomeron loop diagrams potentially 
suppress the enhanced survival probability < I'S'gj^j^l > considerably. However the contribution of higher 
order multiple Pomeron loop diagrams needs to be included, before any precise statements can be made 
about the enhanced survival probability. In the next section, the technique used for deriving the amplitude 
for diagrams with multiple Pomeron loops, will be explained. A good knowledge of how to calculate these 
type of diagrams is a necessary tool for estimating the contribution of the full set of enhanced diagrams, 
to the survival probability. 



4. The main idea 



N = l 



N = 2 





Figure 14: The Pomeron enhanced diagram with 
N = 2 generations of branching. 



Figure 15: The Pomeron enhanced diagram with 
N — 2> generations of branching. 



In this section the approach used to calculate the amplitude of more complicated Pomeron diagrams 
with more generations of Pomeron branching, such as the diagrams shown in Fig. 14 and Fig. will be 
explained. Until now the generally accepted technique for estimating the amplitude for Pomeron enhanced 
diagrams of this type, has been the well known mean field approximation of Mueller-Patel-Salam-Iancu 
(MPSI). In the previous section, the major contribution to the loop amplitude was shown to be equivalent 
to 2 non interacting Pomerons. In ref. [0] we showed that the amplitude of any general Pomeron enhanced 
diagram, is equivalent to the amplitude of the diagram of non interacting Pomerons, with renormalized 
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Pomeron vertices. This was one of the first steps towards deriving the amphtude of more comphcated 
Pomeron enhanced diagrams, other than the simple loop of Fig. |9| in perturbative QCD. The goal of this 
section is to show that the amplitude of an arbitrary Pomeron enhanced diagram, can be estimated in the 
precise QCD approach instead of the MPSI approximation. 



To illustrate some examples of Pomeron enhanced diagrams. Fig. |lj shows the enhanced diagram with 
N = 2 generations of branching, and Fig. ^ shows the enhanced diagram with = 3 generations of 
branching. In Fig. the first generation of splitting forms the larger "embedded" loop, and the second 
generation of splitting forms the two simple loops at the center of the diagram to form a total of 2^ — 1 = 3 
loops. Likewise in Fig. 15 there are 2^ — 1 = 7 loops. Diagrams such as these can be calculated in 
perturbative QCD using the following approach. Consider for example Fig. 16 which shows the diagram 
with N generations of Pomeron branching, leading to a total of 2^ — 1 loops. In this picture the N 
generation diagram is equivalent to 2 sets of diagrams with — 1 generations of branching, embedded in 
one large loop as shown in Fig. Qq. 




Figure 16: Pictorial representation of the Pomeron enhanced diagram with N generations of branching, as the 
product of 2 sets of diagrams each with TV — 1 generations of branching, embedded in one large loop. 

From this picture the amplitude (i(^) (^^1 A, Si/h) of the N generation diagram which occupies a rapidity 
gap A, is a function of the two sets of A^ — 1 generation diagrams embedded in the larger loop, which occupy 
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a rapidity gap Ai and A2 respectively, so that one can write for the amphtude; 



d{N) i^WSyn) = d(^N-i) (I'll^uSyH) ^(tv-i) (i^2|A2) (4.1) 

where denotes all necessary integrations over the conformal variables vi and 1^2, and also over the 
rapidity gaps Ai and A2, where < Ai < A, and Syn < A2 < A. Eq. (|4.lD shows the iterative ex- 
pression which will be used to derive the expression for the arbitrary Pomeron enhanced diagram with 
generations of Pomeron branching. The formalism which will be followed is to start from the basic 
simple loop amplitude {iy\A, dyn) derived in section ^, and plug this into the expansion of Eq. (|4.1|) to 



yield the amplitude d(2) {u\A,SyH) of the N = 2 generation diagram shown in Fig. 14. Similarly inserting 
ci(2) {^\A.,SyH) into Eq. ( |4.1j ) gives the amplitude ci(3) {vlA, Syn) of the N = 3 generation diagram shown 
in Fig. |l^. Continuing with this technique, the general expression (i(7v) A, dyn) for the amplitude of the 
A'" generation diagram shown in Fig. Il^ can be found, using proof by induction. 




Figure 17: The diagram with N generations of Figure 18: (a) shows the contribution to the N 
branching, shown as two sets of iV — 1 generation generation diagram amphtude which stems from the 
diagrams embedded in one larger Pomeron loop. region I: {v, |j/2|} — {0, 1/2, 1/2} and region 11: 

{\v\,vi,V2] = {1/2,0,0}. Regions I and 11 lead 
to the renormalization of the Pomeron intercept, 
and the loops are preserved, (b) shows the contri- 
bution to the amplitude of the N generation dia- 
gram, which is equivalent to the amplitude of 2^ 
non interacting Pomerons, and stems from the re- 
gion{|i.|,|z/i|,|z.2|} = {l/2,l/2,l/2} 



In section ^ the simple Pomeron loop amplitude was calculated for the two regions I: {v, |z^2|} = 
{0, 1/2, 1/2} and II: z^i, z/2} = {1/2, 0, 0}. Region I led to the contribution which is the renormalization 
of the Pomeron intercept (see Eq. (|3.19| )). Region II led to the contribution to the amplitude which 



- 22 - 



is equivalent to the amplitude of 2 non interacting Pomerons, with renormalized Pomeron vertices (see 
Eq. ( 3.26I )). Likewise, calculating the amplitude of the generation diagram shown in Fig. |l^, leads 
to 2 contributions. The first is the contribution from the renormalization of the Pomeron intercept (see 
Fig. 11^). The second contribution is from the diagram which is equivalent to 2^ non interacting Pomerons, 



with renormalized Pomeron vertices, shown in Fig. 18b. These two contributions stem from 3 regions for the 
conformal variables, where in the notation of Fig. 17, these are region I: {v, |i^2|} = {0, 1/2, 1/2}; region 
II: {W\,vi,U2} = {1/2,0,0}; and region III: |i/2|} = {1/2,1/2,1/2}. For the diagrams with iV > 2, 

regions I and II will lead to the contribution shown in Fig. ^a which is equivalent to the renormalization 
of the Pomeron intercept. For N > 2 Region III will yield the contribution to the amplitude which is 
equivalent to 2^ non interacting Pomerons, with renormalized Pomeron vertices. The contributions of all 
three regions to the amplitude of the N generation diagram can be summarized as; 



^(;v) (A, 6yH\Fig. 0) =^L) (A, ^y^lFig. |T| a) + ^{L (A, ^y^lFig. |T8|a) + A]^1 (A, SynlFig. ^h) 



_ a MO) A , ,11 20.(0) A , ,111 2^0.(0) A 

where /c(^) are constants which contain all other terms which are included in the amplitude. 



(4.2) 



5. Multiple-loop Pomeron enhanced diagrams 



In this section, the amplitude for diffractive Higgs production from the Pomeron enhanced diagrams with 
N = 2 and = 3 generations of branching, shown in Fig. |l^ and Fig. |2^ respectively, are derived using 
the iterative technique outlined in section The resulting expressions for these diagrams, will provide the 
tools necessary to form a general expression for the diagram with generations of Pomeron branching. 



shown in Fig. 21 



The scattering amplitude of Fig. IS is given by a straightforward extension of Eq. (|3.12|), namely; 



A(2) (A, ^y^lFig. H) = ^ / di^h {v) (u) e^^'^^^d^^) (^^1 A, Syn) E,E'_,Ah (5.1) 

where d(2) (^ I A, 5yH) labels the contribution to the scattering amplitude of the 2^ — 1 = 3 loops which 
arise from the N = 2 generations of Pomeron branching. The order of the symmetry group of the diagram 



of Fig. |19| is 5'(i)5'(2)) where = 16 is associated with the large outer loop from the same considerations 
discussed in section H], and 5(2) = 8 is associated with permutations of the two internal loops, which lead 
to identical diagrams. 

Following the formalism described in Eq. ([4.1|), d(2) (j^l A, Syu) is a function of the two simple Pomeron 
loop amplitudes d{\) (i/i|Ai, 5y//) and (z^2|Ai) embedded in the large loop as shown in Fig. 19. Here 
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y/ 




QCD color dipole 

BFKL Pomeron 
propagator 

Higgs boson 



Figure 19: The Pomeron enhanced diagram with N = 2 generations of branching, leading to 2^ — 1 = 3 loops. 



Ai = yi — y( is the rapidity gap occupied by the large outer loop in Fig. 19. Hence in this approach, 
including the symmetry factor 1/S'(i)5(2) = 1/2^; 



^ J —00 J —00 

ryi—Syn 

^y/g{^.i) + .M-c.H)Ai^^^^(^^|^^^^y^) d(i)(z.2|Ai) (5.2) 



y 



X / dyi 

'y'+Syn Jy' 



The simple loop amplitude was derived in Eq. ( 3.17| ) and Eq. ( 3.22| ) for two different regions. d{i^ {v2 \ A 



which labels the internal loop in Fig. 19 to the left without the production of the Higgs boson, is derived 
from Eq. ( p. 17 ) or Eq. ( 3.22| ) by setting dyn equal to zero. By inserting Eq. ( |3.17| ) into Eq. (|5.2|), one can 
evaluate the integration over the conformal variables for region I: {v, |z/i|, \v2\} = {0, 1/2, 1/2} and region II: 
{|z/|, z/i, z/2} = {1/2,0,0}. For region I the triple Pomeron vertex derived in Eq. ( |A.22| ) should be inserted, 



and then integrate over over vi and z/2 using the contour C shown in Fig. |13| and sum over the residues 
at {iui,iu2} = 1/2. An extra factor of 2 is included to take into account the identical contribution from 
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the residues at {iui,iu2} = —1/2. Overall after integrating over the rapidity variables, the contribution of 
region I to Eq. ( ^ ) is given by the expression; 



4) (HA, dyn) =d'<f Q + iv^ ' Q - X (^) e-^^^^'y^ 



(5.3) 



where <i' = ^ 1 1 ) (54) 



Inserting Eq. ( |5.3| ) into Eq. ( |5.1|) and integrating over using the method of steepest descents in the 
same way described in section ^ yields the following contribution of region I to the scattering amplitude 
for Fig. Ill 



^^2) (A,(5y^^|Fig. [1) 



„a.(0)(A-5s//f) 



1287r2iV2 



XX (0) 



(.."(0) {^-6yH)f'^ 

u? (0) (A - byu) 



6 - -^a^byu ) (0) ( w (0) + 



(5.5) 



8.91 X 10"^^ GeV-2 
3.3 X 10"^^ GeV"2 



{a, = 0.12) 
(os = 0.2) 



(5.6) 



Alternatively, after inserting Eq. ( p. 17 ) into Eq. ( |5.2| ), the contribution from region II, namely i^i, z/2} 



{1/2,0,0} is found by substituting for the triple Pomeron vertex Eq. (A. 23) and integrating over the con- 
formal variables vi and U2 using the steepest descents method, which yields the expression; 



a'd 



I j2 



{l/2 + iv){l/2-iv) 



1 



dyi 

y'+Syn •> y 



xx' {^)\oJ (0) + -vo' (0)(Ai - dyu) (0)+ -uj' (0) Ai 



yi-Syn 

dyi 
1 



,(2w(0)-cj(i/))Ai 



(Ai (Ai - 6yH)) 



3/2 



where a ' 



az 



267riV4[^//(0)]3 



(5.7) 
(5.8) 



Inserting Eq. (5/7) into Eq. ( |5.l[) and using the contour C shown in Fig. 13 for the z/ integral, the 
solution is the sum over residues at iu = 1/2. In the same way an additional factor of 2 includes the 
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identical contribution from the sum over residues at ii^ = —1/2. Using this approach, the contribution of 
region II to the scattering amphtude of Fig. |l^ is; 



A" ( A, 6yH\Fig. m = (0) (A - Syn) 



a, dA 



297riV2 



4.94 X 10 



-17 



GeV" 



3.23 X 10"^^ GeV 



T-2 



{as = 0.12) 
(a, = 0.2) 



(5.10) 



Eq. (5^) and Eq. ( |5.9D are the contributions to the scattering amphtude of Fig. 1£, which lead to the 
renormalization of the Pomeron intercept. Alternatively, the contribution to Fig. |l^ which is equivalent to 
4 non interacting Pomerons, is found by substituting for ci(i) (z^i| Ai, (^y//) and (z^2|^i) the simple loop 
amplitude of Eq. (3.22) in Eq. (5.2), and evaluating the integrals over the conformal variables, taking into 
account the singularities which arise from region III: {|z^|, |i^2|} = {1/2, 1/2, 1/2}. Inserting the triple 



Pomeron vertex of Eq. (A. 22) into Eq. ( [5. 2] ), the ui and U2 integrations are evaluated by summing over the 
identical residues at {iui,iv2} = 1/2 and —1/2, which arise from integrating along the contours C and C 



shown in Fig. 13, yielding the result; 



d^^{u\A,6yH) =aHe-^^^^'y" {1/2 + iuf {1/2- iufx{^) dyi dy/e-^^'^)^^ 

Jy'+SyH Jy' 

x(Ai-5y^,)Ai — -— — — — ^ U —r^\ (5-11) 




as dAi 



A? 



where b 



at 



1 



1 



(5.12) 



210 Ar2^ 

Inserting Eq. ( |5.1l| ) into Eq. ( |5.l[) , the u integration is solved using the same method of residues 



described above, so that the contribution of region III to the scattering amplitude of Fig. |19|, which is 
equivalent to the amplitude of non interacting Pomerons is; 



A^^{A,5yH\Fig.^ 



X (A - Synf A 



2'^N^TT 

7-1 dV eM0)A 
\asdA) A^/'^{A-5yHf'^ 




2.85 X 10"^^ GeV-2 



1.39 X lO"'^ GeV 



T-2 



[as 
{as 



0.12) 
0.2) 



-1 d y e^o)^ 

^dA j A3 



(5.13) 



(5.14) 
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Prom an observation of Eq. ( 5.13 ); 



^™ {A,SyH\Fig. m 



\2) 



(5.15) 



where ^(2) contains all the other terms contained in the amplitude. The above results show that 



the main contribution to the the scattering amplitude of Fig. 19 is just the equivalent of 4 non interacting 
Pomerons, with renormalized Pomeron vertices. The compete scattering amplitude of the N = 2 generation 
diagram of Fig. 19, is the sum over the contributions of Eq. (^.5|), Eq. (|5.9|) and Eq. (5.13D, namely; 



A{A,6yH\Fig. H) = A) iA,6yH\Fig. li + (A, <5yj^|Fig. H) + A™ {A,6yH\Fig. |D (5.16) 



: 2.85 X 10"^^ GeV-2 (a, = 0.12) 
1.39 X lO"'^ GeV-2 (a^ = 0.2) 



(5.17) 



Comparing the values of Eq. ( 3.2^ ) and Eq. ( |5.17| ), Fig. 1£ is of the same order of magnitude as the 
simple loop amplitude of Fig. |^. This implies that in order to estimate the enhanced survival probability, 
it is not enough to just take into account the simple Pomeron loop, but rather also enhanced diagrams 
with multiple Pomeron loops need to be included in the estimate. 



Fig. 20 shows the diagram for diffractive Higgs production in t-channel Pomeron exchange, where there 
are 3 generations of Pomeron branching, which recombine to form 1 large loop at the = 1 level, 2 smaller 
internal loops at the N = 2 level and 4 simple loops at the N = 3 level at the center of the diagram. The 



scattering amplitude of Fig. 2C is given by the expression; 



^(3)(A,5yH|Fig.|g) =-^/ d,^hiu)X^iu)e'-(-^^d^^){iy\A,5yH)E,E'_,AH (5.18) 

where d(3) {A,6yH) is the contribution of the 2^ — 1 = 7 loops in Fig. ^to the scattering amplitude. 
Following the approach of Eq. (|4.l|) , d(3) (A, Syn) is a function of d(2) (Ai, dyn) and d(2) {Ai,5yH) embedded 
in the larger outer loop, as shown in Fig. |2^. The = 3 generation amplitude d(3) (A, dyn) should also 
have in the denominator the same symmetry factor 1/2^ as for the case of the N = 2 generation amplitude, 
for the same above explained reasons. As such d(3) (A, Syn) is given by the expression; 
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di'ih (I'l) X"^ (i^i) / dv2h {^2) X^ {v2) 



dyi 
y'+^VH Jy 



dy;e(^M+-M-^M)^i d(2) ( I AuSyn) d(2) (^^2 I Ai ) (5.19) 



-27- 



Figure 20: The N—3 generation Pomeron loop diagram, leading to 2^ — 1 = 7 loops. 



Following the same steps as Eq. (|5.2| ) - Eq. ( 5.13| ), the contributions of region I and II to Fig. ^ 
which leads to the renormalization of the Pomeron intercept, and the contribution of region III which is 
equivalent to the amplitude of non interacting Pomerons are given by the expressions; 



^^3) (A,<52/^^|Fig. |0D 



1287r2iV2 
x6 ( 6 - ^as^VH 



16 4 



x(o) 



9 



{2as)SyH]oj' (0) L^(0) + 



UJ 



{0){A -6yH) \ 



(5.20) 



3.93 X 10"^^ GeV"^ 
6.27 X 10"^^ GeV'2 



{as = 0.12) 
{as = 0.2) 



(5.21) 
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AjJ (A, 5yH\Fig. m= "^"X'^^'^ e^e - lasSyn'^ (0) co' (0) (A - Syn) 



-1 dV \ e2^(Q)(^-^?/g/2) 

\(A(A-5yH))^/^ ^J^U + 

= 2.52 X 10^26 GeV~2 
8.53 X 10^2^ GeV"2 



(0)(Ai-(^j/j, 
3 

[as = 0.12) 
(as = 0.2) 



(5.22) 



(5.23) 



A™ (A,(5yH|Fig. 1^) 



29A^2^ 



xA^ {^-5yHY 



-1 d Y 



„2a.)(0)A 



a. 



3.34 X 10-2^ GeV"2 
1.28 X lO"'^ 



dAy A3/2(A-5y^ 

(as = 0.12) 
{as = 0.2) 



,3/2 



a, dA 



A3 



(5.24) 



(5.25) 



Eq. (5.24) is the contribution which is equivalent to the amphtude of 8 non interacting Pomerons, as 
can be seen clearly by recasting Eq. ( ^.24 ) in the following form; 



^™ (A,<52/H|Fig.|2C 



(5.26) 



where /c(3) contains all other terms included in the amplitude. The major contribution comes from the 
non interacting Pomeron contribution of Eq. (|5.25| ). The full scattering amplitude of Fig. 20 is the sum 
over the contributions of Eq. ( |5.20| ), Eq. ( 5.22| ) and Eq. (5.24), namely; 



A (A, 5yH\Vig. ^ = ^Jj) (A, 5yH\Vig. ^ + a}1 (A, 5y/,|Fig. ^ + A}]] (A, SyuWig. 0) (5.27) 



H3) 



\3) 



: 3.34 X 10^2^ GeV-2 (a^ = 0.12) 
1.28 X lO"'^ GeV-2 (a^ = 0.2) 



(5.28) 



The value in Eq. ( 5.28| ) is very close to the amplitude of Fig. ^ for = 0.2 (see Eq. ( |5.17| )) and the same 
order of magnitude as the amplitude of the simple Pomeron loop of Fig. |^ (see Eq. ( |3.28D ). This shows 
that in order to obtain an accurate estimate of the enhanced survival probability < | S'^^Yi I > i sum over 
the complete set of Pomeron loops needs to be found. This is the subject matter of the next section. 



- 29 - 



6. The summation over Pomeron loop diagrams 



In this section the diffractive Higgs production amphtude from the general multiple Pomeron loop diagram, 
with generations of branching is derived. This general expression provides the basis for summing over 
the complete set of Pomeron loop diagrams, in perturbative QCD, instead of using the MPSI approach. 
In a diagram with generations of Pomeron branching shown in Fig. the A;*^ generation of branching 
gives rise to l:^"^ loops, so that in Fig. 21, there are a total of ^Y^=\'^~^ — — 1 loops. So far 3 
contributing regions to the amplitude of Pomeron loops have been discussed. An analytical formula for 
the diagram with A^ generations of branching can be derived, for all 3 of these regions using the iterative 
technique outlined in section ^, without any a priori assumptions. For the diagram of diffractive Higgs 
production in t-channel Pomeron exchange, where there are A^ generations of Pomeron branching as shown 



in Fig. 21, the high energy scattering amplitude is labeled; 




Figure 21: The diagram with N generations of Pomeron branching, leading to 2 — 1 loops. 



yl(^)(A,<5yH|Fig.il|) =^ / dvh{v)\^(y)&'^''^^di^^^{y\AMH)E,E'_,AH (6.1) 



where the amplitude for the A^ generations of loops d(^) ('^l^i ^Uh) is found from the iterative equation; 



-j^ f'QO poo 

^ -'— oo J —oo 



Jy'+Syu Jy' 

xd^N-i)M^i,5yH) d^N-i){^2\^i) yN>l. (6.2) 

The expression for A, Syn)-, d(2) (^^1 A, Syu) and d(3) Syn) for the diagrams of Fig. ^, Fig. ^ 

and Fig. |2^, were derived using the iterative expression of Eq. (^). Continuing this process using identical 
techniques, from observing the outcome expression for ci(jv) (j^l A, dyn) for higher values of N , a predictable 
sequence emerges. From an observation of this sequence, a general formula for (i(jv) (^^1 A, can be 
written for the contribution of region I: ( {u, |f2|} = {0, 1/2, 1/2} ), region II: 

( i/i, 1^2} = {1/2, 0, 0} ) and region III: ( |i/2|} = {1/2, 1/2, 1/2} ). This expression can be 

proved by induction, whereby provided the formula for d(jv) (i^|A,(5yH) is true for N = 2, and yields 
the predicted formula for + 1 by plugging d(7v) (z/|A,(5yH) into Eq. (|6^), then the proof is complete. 
In this formalism the following expressions for the contribution of the above three regions are derived for 
d{N) i^l^^Syn) ; 



^ (k-2)2l^-^] f k + I (A; + 3)(A: + 2)(A: + 1)^ 



X a 



2k^ 



n 

k=2 

T^/fc + l (A; + 3) (A: + 2) (A; + 1) A; + 1 ^ 
11 1 7^ ^ j:^asOyH 



k=2 



2A;2 



A;2 



VA^ > 1; (6.3) 



("l^-*") = (l/2 + i.) (l/2-i./ '^ - l-^^)^*^ - 1' 

/■J/ ryi-Syn g(2w(0)-a;{i/))Ai 

X / dyi / dy/ —pr 

Jy'+Syn Jy' (Ai (Ai - y//)) ^ 

xx=(0).--*(0)(.(0) + ii!(21<^P«)(w(0) + ^:^) ViV>2; (6 A) 
where the function / (A^, dyn) is defined by; 
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f{N,6yH) = W a 

k=2 
N 

k=2 



(fc-2)2[^-^l (k + l _ (fc + 3) (fc + 2) {k + ir ^'""""-^ 
k 2A;2 



n 



k + l (fc + 3) (A; + 2) (A; + 1) + 1 _ ^ , 
— 2k^ + ; 



(6.5) 



and / {N — 1) = / (A^ — 1, Syu = 0) can be read off Eq. (S.5) by setting Syu = 0, and 



{ah) 



2[JV-i] 



-1 d 



a.. 



y'+Syu Jy' 
p2cj(0)Ai 



'"X'Af-^-^ (Ai - Jy^)^-^ e-M^^ 



dAi; A?/'(Ai-5yj^)3/2j\V«. dAi; A? 



1 d \3eM0)Ai- 



2['v-il-l 



ViV > 2. (6.6) 



wliere tlie full set of constants are; 



29a! 



2 



A^>[u;" (0)]3' 



d' 
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aZ 
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26iV47r[w"(0)]3^ 



210 V m 



(6.7) 



Inserting the formulae of Eq. ( |6.3| ), Eq. ( |6.4| ) and Eq. (6.6) into Eq. ( |6.1D , one finds the following 
expressions for the contributions of regions I, II and III to the scattering amplitude for diffractive Higgs 



production, from the N generation diagram of Fig. 21 
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ViV > 2; (6.10) 



Eq. ( |6.8D and Eq. ( |6.9D are the contributions to the amphtude of the N generation diagram of Fig. ^ 
which leads to the renormalization of the Pomeron intercept. Eq. ( 6.10 ) is the contribution to Fig. ^ 
which is equivalent to the amplitude of 2^ non interacting Pomerons, with renormahzed Pomeron vertices. 
This can be seen by observing that Eq. ( |6.10| ) can be written in the form; 



^2^a;(0)A 



(6.11) 



where A:(jv) includes all other terms contained in the amplitude. The complete expression for the 
scattering amplitude of Fig. ^ is the sum of the contributions of Eq. ( |6.8D , Eq. ( |6.S| ) and Eq. ( |6.10 ), 
namely; 
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^(;v) {A,6yH\Fig. H) 

= ^J^) {A,6yH\Eq. (^) + A^},^ {A,6yH\Eq. Q) +^J]^) (A,5yH|Eq. {^) (6.12) 

Table || lists the results for the multiple Pomeron 
loop amplitudes for = 0.12 and = 0.2. The 
greatest contribution to Pomeron enhanced diagrams 
stems from region III, which leads to the contribution 
to the amplitude which is equivalent to the amplitude 
of non interacting Pomerons. The values in the table 
indicate that the amplitude of the N generation dia- 
gram becomes smaller as grows, for energies within 
the LHC range. The difference in values for different 
values of a<j show that the amplitude is very sensitive 
to the choice of the Pomeron intercept. The sum over 
Table 1: Results for the scattering amplitude for Pomeron loop diagrams from iV = (so that the ba- 
diffractive Higgs production from the multi Pomeron sic diagram of Fig. || is included), takes the following 
loop diagram with N generations of Pomeron branch- form; 
ing, for the contribution of regions I , 11 and III. The 
mass of the Higgs boson is assumed to be Mh = 
100 GeV, and the rapidity gap A between the scat- 
tering protons is taken to be A = 19, based on proton 
proton collisions at the typical LHC energy y/s = 14 
TeV. 





A^N) (A = 19, 6yH = In (M^/4so)) 
Us = 0.12 as = 0.2 


iV = 


1.57 X 10"^ 2.17 X lO"'^ 


N = 1 


2.05 X 10-1° 2.15 X lO-'^ 


N = 2 


2.85 X 10-^^ 1.39 X 10"^ 


N = 3 


3.40 X 10^24 28 X 10"^ 


N = 4 


8.20 X 10-^2 2.31 X 10-*^ 


N = 5 


8.13 X 10""^^ 5.11 X 10'^ 


N = 6 


1.36 X 10-1^6 4.26 X 10-1° 



^{A,6yH) =^(o)(A,(^yH|Fig.0; Eq. (p|)) + ^ {-if A^^-j {A,6yH\Eq. (^)) (6.13) 



N=l 



Fig. 22 shows the plot for the amplitude of the single Pomeron diagram derived in Eq. ( p. 24 ), 



for the basic diagram shown in Fig. |^ (upper line), next to the plot for the sum over Pomeron loops 
S (A,(5yH|Eq. ( |5.13[ )) up to = 20 (lower line), against the rapidity separation A between the incoming 
projectiles. From an observation of Fig. ^ it is clear that for rapidity values approaching the typical 
LHC range A = 19, S (A, 5yj|/|Eq. ( |6.13D ) starts to be substantially less than the basic single Pomeron 
amplitude, as the two graphs grow further apart. This proves that the shadowing correction of Pomeron 
loops to the basic diagram of Fig. ^, is large within the LHC range of energies. As the rapidity separation 
between the scattering protons increases, the effect of this shadowing correction increases. 
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Figure 22: Plot for the single Pomeron amplitude A(o) (A, SynlFig. ^ derived in Eq. (2.24) (upper line), next to the 
plot for the sum over Pomeron loops E (A, Synl Eq. ( 6.13| )) up to = 20 (lower line), against the rapidity separation 
^ = y — y' between the scattering protons. The values of A go from Syn = In (Af^/4so) (assuming the Higgs mass 
Mh = 100 GeV and sq = 1 GeV^), up to the typical LHC rapidity A = 19. as = 0.2. 



7. The survival probability for exclusive Higgs production 



This section is limited to an estimate of the contribution of enhanced diagrams to the survival probabihty 
< l^g^jjl > of large rapidity gaps, in exclusive diffractive Higgs production shown in Fig. ^. In section ^ the 
sum over Pomeron enhanced diagrams in Eq. ( |6.13| ) was derived in the exact QCD approach, thanks to the 
expressions of Eq. ( |6.8D , Eq. ( |6.9| ) and Eq. ( |6.10| ) . This means that the complete set of hard re-scattering 
contributions to the enhanced survival probability can be estimated in QCD to any order, without relying 
on the mean field approximation approach. 



The enhanced survival probability < l^g^j^l > is estimated by subtracting from the basic diagram of 
Fig. ^, the first enhanced diagram of Fig. ^, and subtract from this the second enhanced diagram of Fig. H 
and so on. Finally divide by the amplitude of the basic diagram of Fig. ^ to give the correctly normalized 
survival probability as the expression; 



< \S. 



cnh 
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A(o) (A,5yj^|Fig. 0) - (A,5yj^|Fig. |) + ^(2) (A,5yj^|Fig. 0) 



A,m (A,(5y^^|Fig.|0D 



N 



N=0 



^(0) {^,5yH) 



(7.1) 
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where yl(o) (^) ^UnlFig. 0) was calculated in Eq. ( ^.24 ) and all the subsequent terms ^(at) (A, dyn) y N > 
1 are the sum over Pomeron loops given in Eq. ( |6.8D - Eq. ( 6.12| ) 




Figure 23: The enhanced survival probability < > plotted against the rapidity separation A of the scattering 

protons, as = 0.2. 

The results for the enhanced survival probability < I'^g^j^l > are shown in the graph of Fig. 23 against 
the rapidity separation A between the incoming projectiles, where all the terms in Eq. (^^) up to = 20 
have been taken into account. The results of Fig. ^ show that < ISg^^j^l > is small, and could be even less 
than 1% for the LHC range of rapidity values, in agreement with refs. The observation from Fig. ^ 

that < ISgjjjjl > decreases as the rapidity gap A between the scattering protons increases, also matches with 
the findings of refs. ||3|, ^. However Fig. 23 shows a steeper rise in < 15*^^^^^! > as the rapidity separation 
A decreases than the results found in refs. ^, where the improved Mueller-Patel-Salam-Iancu (MPSI) 
approach was used. We believe this difference in the slope of the graph for < l^g^^j^l >, shows that there are 
more complicated diagrams to take into account other than those diagrams considered here, which would 
lead to a result for < IS'^j^j^l > which shows a more gradual increase as the rapidity gap A decreases. This 
will be discussed in more detail in the conclusion section. 



8. Conclusion 

The main achievements of this paper are the following; 

1. The observation that the major contribution to the Pomeron loop amplitude stems from the non 
planar piece of the Korchesky triple Pomeron vertex of Fig. ^ 

2. The derivation of the general analytical expression for the diagram with an arbitrary number of 
Pomeron loops, found in Eq. ( |6. - Eq. ( |6.12| ). 

3. The finding that the main contribution to the amplitude of multiple Pomeron loop diagrams, is 
equivalent to the amplitude of non interacting Pomerons. 
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4. The solution to the summation over Pomeron loop diagrams, for the first time in the exact QCD 
approach. 



5. Estimation of the contribution of Pomeron enhanced diagrams to the survival probability < |5'g^j^| > 
for diffractive Higgs production in the QCD formalism, instead of the MFA approach. 



The main achievements of this paper are the formulae of Eq. ( |6.§| ) - Eq. ( 6.12| ) for the amplitude of 



the multiple Pomeron loop diagram, with A'^ generations of Pomeron branching as shown Fig. 21. It was 
found that there are 3 regions which contribute to Pomeron loop diagrams. Regions I and II lead to the 
renormalization of the Pomeron intercept. Region III leads to the contribution to Pomeron loop diagrams, 
which is equivalent to the amplitude of non interacting Pomerons, with renormalized Pomeron vertices. The 
contribution from non interacting Pomerons, gives the greatest contribution to the amplitude of Pomeron 
loop diagrams. The property of the equivalence of Pomeron loop diagrams to non interacting Pomerons, 
was first noticed by A. Mueller (see ref. |5|). We showed in ref. 0, that if the sum over Pomeron loops 
was performed in the QCD approach, then the summation should reduce to the sum over non interacting 
Pomerons diagrams. Therefore the findings in this paper are in agreement with the results in the above 
mentioned papers. 



The sum over Pomeron loop diagrams yields a significant shadowing correction to the basic single 
Pomeron diagram of Fig. ^. The shadowing correction however is only significant for a rapidity separation 
between the incoming protons A > 14, and the shadowing correction becomes larger as the rapidity sepa- 
ration grows (see Fig. |2^). We did not take into account the diagrams which contribute to the vertex in 
the framework of the Schwinger-Dyson equation. However, since the intercept of the Pomeron A > as 
was shown in refs. ||5|, 0, these diagrams give a small contribution. 



The application of the sum over Pomeron enhanced diagrams to an estimate of the enhanced survival 
probability < {S'^^^il >i produced a result which shows that < l^g^^j^l > is potentially less than 1% for a 
typical LHC rapidity separation A = 19. From an observation of Fig. < l^g^^j^l > decreases as the 
rapidity separation A increases (see Fig. 23). Both of these findings are in agreement with refs. ||3|, ^. This 
comparison shows that the summation over Pomeron enhanced diagrams in the formalism of this paper, 
reproduces the expected behavior of < IS^^j^J > with rapidity separation A. The estimates for the LHC 
range of rapidity values are close to the previous estimates of refs. M, 0] (i.e. less than 1%). 



The behavior of the enhanced survival probability as a function of the rapidity separation, is a property 
which stems from the definition of the survival probability. That is the survival probability is a quantitative 
measure of the effect of shadowing corrections, that stem from inelastic scattering emerging from extra 
parton showers in the reaction proton + proton — > proton + [LRG] + Higgs + [LRG] + proton, (where 
[LRG] denotes a large rapidity gap). Increasing the energy of this reaction leads to a rise in the number 
of extra parton showers, which spoil the large rapidity gaps as shown in Fig. ^ This leads to the natural 
conclusion that as the rapidity separation between the scattering protons increases, the survival probability 
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decreases as more unwanted parton showers arise. 




WAmW BFKLPomeron 

ly^y^y^^ Renormalizcd BFKL Pomeron 

, Higgs bofion 



Figure 24: More complicated Pomeron enhanced Figure 25: Schematic representation of loop correc- 
diagrams that have not yet been taken into account. tions that leads to the renormalized BFKL Pomeron 
(a) shows a diagram with N — 2 generations of propagator. 
Pomeron branching, (b) shows = 3 and (c) shows 
N — 4: generations of branching. 



The slope of < |5'gjjjj| > of Fig. |^ with the rapidity separation is steeper than the findings of refs. 
IP, Q], which relied on an improved (MPSI) model. This observation we believe shows that the summa- 
tion over Pomeron loops in this paper, while this is an important step forward, is possibly not yet complete. 



For example, more complicated enhanced diagrams such as those shown in Fig. 24 could provide an 
important contribution to the enhanced survival probability, which would lead to a behavior with a less 
steep increase with rapidity separation. The diagrams in Fig. ^ are not topologically equivalent to any 
of the other enhanced diagrams discussed in this paper. Since multiple Pomeron loop diagrams decrease 
with the number of generations of branching N, diagrams (a), (b) and (c) in Fig. ^ which belong to the 
N = 2, N = 3 and = 4 class of loop diagrams are not expected to be larger than the simple Pomeron 
loop diagram of Fig. S, however their contribution to the summation over enhanced diagrams could affect 



the slope of Fig. 23 to make it less steep. There is also the issue of the renormalization of the BFKL 



Pomeron which needs to be taken into account. For example diagrams such as Fig. |25|, which although are 
not expected to be large for any value of could still provide an important contribution to < l^g^^j^l >, 
which would make the slope with rapidity separation not as steep. 

We hope that the findings of this paper will be useful for experiments aimed at discovering the Higgs 
boson at the LHC. We also hope that the methods developed in this paper will provide an important 
foundation for future calculations of BFKL Pomeron loop diagrams. 
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A. The triple Pomeron vertex 

In this section of the appendix the derivation of the triple Pomeron vertex, which was first done by 
Korchemsky in ref. M will be outlined. The triple Pomeron vertex is the vertex of three BFKL Pomeron 



states, which couple either as shown in the planar diagram of Fig. 11, or they couple in the way shown in 



the non planar diagram of Fig. 12, The coupling of three BFKL Pomeron states with the center of mass 
coordinates R, Ri and R2 and the conformal weights 7, 71 and 72 is given by for the planar and non-planar 
diagrams respectively; 

n {y I yi,v2\Fig. O) = / j d^^i J "^^^^ \Rq^\^\iI^\^\R^q\2 (^'^^ 
A {u I i^i,U2\Fig. |12|) = y d^Ro j d^^ijj^ ^7 ^71 ^72 (A.2) 

The leading contribution to the high energy scattering amplitude stems from the region where all the 
conformal spins vanish, namely n = ni = 77-2 = 0. Adding together Eq. ( [A.l| ) and Eq. ( [A. 2 ) and integrating 
over the complex coordinates, the triple Pomeron vertex takes the form; 



r(Hz.i,i.2) = {4z.2 + i}2 (^f^(^,|^,^,^2) + |^AH^^i,z^2){xM-x(z^i)-x(^^2)}) (A.3) 

where x (^) is the function; 



X{u) = ^e{^l^{l)-i^{l- + iu]} (A.4) 



2 



and where 



n{u\u^,u2) = vr^ jr^Q + .z.) r^Q + iz.!^ ^{l''") ^{\~'''^ 

X J'' ( 1/ , z^i , 1/2 ) Ja ( , z^i , 1/2 ) (A.5) 
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where the contracted a index implies summation over a = 1,2,3. The integrals Ja{v , i^i , 1^2) and 
Ja {i^ , 1^1 , 1^2) were evaluated in ref. where they were found to be; 



Ji{u,ui,U2) = T + iu + ii/i - iu2^ T ^ ^ + ii/^ T ^ ^ - ii/^ T ^ ^ + ii^i ^ F ^ ^ - ii^i 
X j dx {1 - x)~^^'^~''''^ 2F1 iu^llx^ 2F1 (^+11^1,^- iui;l\x^ 

J2 {u, ui, V2) = T + iu2 ^ ^ Q + ^'^^ ^ ~ ^ ^ ~ 

(i - ^.2) 



T {\ + iv — 11/2) T (1 — iu — ii'2) 

/I 1 1 1 , 

X4,Fs\-+il^l, --ivi, --tU2, --IV2\1 + IV-IV2, I - IV - IV2 , 1 | 1 



Jl {V, Vl,V2) 



r (1 + ivi — iu2) r (1 — ii^i — ii^2) 

^ A 1 1 1 

,F,[^- + ^u,--^u,--^.2,-- 



/I 1 1 1 , 

X A^F^ \-+ IV , -- IV , -- IU2 , -- IV2\1 + IVl - 11^2 A - - ^^2 



hiv, Vi,V2) 



J3 {V, Vl,V2) 



r (i - iv2) r (1/2 -iv- ivi + iv2) 

X J\xx (1 - ^)-l/2+-2-..-..i Q - " ^1"^) ''''' ^'^) 

r (1 + iu,) r (1 + r - ii^) (i - iz^i) 

r (1/2 — zi/ — ivi + ii/2) 

X jjl' (1 - Q _ 1 _ ^p^ Q _ 1 _ 

r(i+i^)r(i + ir.2)r(i-zz.i)r^(i-iz.) 

r (1/2 — iz^ — ivi + iz^2) 

X (1 - ^)-i/2— — Q - ^ - ii^; l|x) 2i^i Q - ii^i, I - ivi; l|x) (A.6) 



For the non planar piece; 
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2 r (i - iz/2) T + + iui - 11/2) 



r + ii/) r (I + zi/i) r (1 + — iz/2) t {1 — iu — 1^2) 

Mil, 

- + iz^i , - - iz/i , - - ii'2\l + iv - iiy2 , I - iv - 11^2] 
k z z z 

r (i - iz/i) r (i - ziv) r (i + + iui + zz^2) r - + ii/i + iv2) 



T {l-iv + iz/2) r (1 + ivi + iz/2) 
X 3-^2 ( ^ + , ^ ~ ' ^ - + ^i^i + «i^2 I 1 - + ^^^2 , 1 + ii^i + i;/2 I 1 ) (A. 7) 



In ref. Q it was shown that ^{i>\i>i , 1/2 ) and A (1/ 1 i^i , 1/2 ) are symmetric under cyclic permutations 
of {u , h'l , 1^2}- That is; 



, U2) = Q{l'2\l' , i^i) = il, {ui \ V2 , u ) 

A (zv I z^i , z^a ) = A (z/2 I z^ , z^i ) = A (i/i | zy2 , ) (A.8) 



with an identical result for the complex conjugates. For specific values of the conformal variables 
{i^, 1^1,^2,}, inserting the appropriate values into Eq. ([A.6|); 



1 1\ ^ 0^3 (1 - ^^1 - ^^2) 

'2'2y' (l/2-ii.if(l/2-iz.2r 



n{iiy\iui,iu2) = n{iu\-,-] = 2ir^ ^ 2 (^-9) 



«>(HM,i.,) = f!(H-l-l) - (l/2 + J)tl/2-..) ^K'''''''-'''G^''^)) 

n(ii/iui,iu2) = n{ -0,0] =— + ^3^2 <^ - , - , - 1, - 1 ^ A.ll) 

1 \ 



Similarly from Eq. ( |A.7| ); 
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2 2y (1/2 - zi/i) (1/2 - zi^2j 

A*(ii/|izyi,izy2) = A(ii^|-^,-^) =7777— ^T7Ta3 ^ ^^-^^^ 

11) [1/ 1 + iv) (1/ 1 — lU) 



1, \ 7r2 

-0,0 = 

2' y (1/2 -ii/) 



A(ii/|zzyi,izy2) = A ( -|0,0 ) =771:^ —n (A.15) 



1,„ „\ vr2 



A*(Hii'i,i^2) = A(--|0,0 = . .2 (A.16) 

1 J [1/2 + 



Collecting the above results and inserting them into Eq. ( [A.3| ), one finds that the planar piece gives 
a first order pole at {v, \ui\, |i^2|} = {0, 1/2, 1/2} and the non planar piece is finite in this region, so that 
the planar piece is the dominant term for this region. For the region i^i, 1/2} = {1/2,0,0} one finds 
the reverse, namely the planar piece is finite whereas the non planar piece offers a first order pole at 
{|z/|, z/i, z/2} = {1/2, 0, 0}, so that the non planar piece is the relevant part of the triple Pomeron vertex for 
this region; 



r {ii'\ivi,iv2) = r I ii^ 



1 1 

'2' 2 



r* {iv\ivi, iu2) = r yiv 
T {iu\iui,iu2) = r f^|0, 



r* {iu\iui,iu2) = r ( — -|0, 



1 1 
2'~2 



61 TT ttg ,^ ,„ 



(1/2-ii/i) (\ll-iv2) 
l-iui-W2-2/Nl 2 
(1/2-ii/i) (1/2-^1/2) iV^^^i^J 



+ iu 



1 



1/2 -iu 
N'i 1/2 + iv 



(A.17) 

(A.18) 
(A.19) 

(A.20) 



The following notation is used for the product of two expressions of the triple Pomeron vertex with 
specific values, where one is the complex conjugate of the other; 



T {iv\ivi,iv2) X V* {iu\ivi,iv2) 



r (z/|z/i, U2) 



(A.21) 



In this notation, combining Eq. ( A. 17] ) and Eq. ( |A.l^ ) yields; 
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32 ( 2TTf at I 1 
1 — ivi — iu2 



1 

nI 

2/iVe' 



{1/2 -iui) {1/2 -iu2) 



{1/2 -iUi) {\/2-iV2) iV2 

and collecting the results of Eq. ( [A. 191 ) and Eq. ( [A.2[l| ) ; 



(A.22) 



r {v\v\, U2) 



B. The fan diagram 



10,0 




^2TT^al 



WWWv 



QCD color dipole 
BFKL Pomeron 



Higgs boson 



Figure 26: The first fan diagram with N 
tions of branching. 



1 genera- 



[l/2 + iv){l/2-iv] 



+ finite terms 



(A.23) 



In this section of the appendix the amplitude of 



the diagram of Fig. 26 known as the Pomeron "fan" 



diagram is calculated. The first fan diagram is shown 



in Fig. 26, where the t-channel Pomeron splits into 



two daughter Pomerons, and the Higgs boson is pro- 
duced from one branch of the fan diagram. The scat- 
tering is between two QCD dipoles separated by a 
rapidity gap A = y — y' in rapidity space, and the 
Pomeron vertex where the splitting occurs is at ra- 
pidity yi where y' + 6yH < Vi < y- The motiva- 
tion for this inequality condition is to ensure that the 
scattering is energetic enough to produce the Higgs 
boson, and therefore the separation in rapidity space 
between the vertex at the top of the fan and the target 
dipole should be at least Syn- In the representation 
of complex angular momentum oj = 'i + j the fan 



diagram amplitude of Fig. 26 has the expression; 



F {u,ui\Fig. li = i (y)'^^^ / / ^^2 5^ (7)50.1 (7i)5^-c.i (72) E^E^.E^^ 

X j cPR j (fRi J d^R2V {R,Ri,R2h ^ ^1,^2) (B.l) 

where following all the conventions of section |2|, the (7) are the propagators of the three BFKL 
Pomeron states in Fig. ^ and the E^ represent the couplings of the BFKL Pomeron propagators to the 
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QCD color dipoles. S is the order of the symmetry group of Fig. Including the 2 ways of attaching 
the Pomeron branch lines at the vertex, and the 2 permutations of the reggeized gluon lines in the triple 



Pomeron vertex shown in Fig. IC , the total symmetry factor in Eq. ( |B.1| ) is S" = 4. The reason behind the 
coefficient (as/4)^ in Eq. ([B.lD follows from similar arguments given in section 0. In Fig. 26 there are three 



vertices where the BFKL Pomeron couples to the QCD color dipole, namely at the top of the diagram 
and twice at the lower ends of the two fan branches. Each vertex brings a factor of ag on account of the 
2 couplings of each of the 2 reggeized gluons in the BFKL Pomeron structure (see Fig. |5|). Each Pomeron 
- dipole vertex contains a 4 - fold degeneracy, since both of the 2 reggeized gluons can couple to either of 
the 2 quark lines in the color dipole, making a total of 4 identical possibilities for each vertex. Thus one 
arrives at a factor of as/4 for each of the three vertices in Fig. |2^, which ensures that identical diagrams 
are not counted more than once. 

y (i?, -R2I7 — *■ 7i;72) denotes the triple Pomeron vertex for the splitting of the BFKL Pomeron 
state labeled by the curly brackets (7 is the conformal variable and R is the center of mass coordinate) 
{i?,7}, into the two BFKL Pomeron states labeled by {i?i,7i} and {i?2;72} (see Fig. |l^). Recall from 
section ^ , that from the condition of conformal invariance comes from the property for the triple Pomeron 
vertex UlTl 



ViR,Ri,R2h - 71,72) = Ro^""'' R^2'''' Ki"^'' R*i2^'' R*20^'°'^ ^ ^2) (B.2) 

where Roi = R — Ri {i = 1,2) and R12 = Ri — R2, and where for example Aqi = 7 + 71—72; and 
^12 = 71+72 — 7) with an equivalent definition for the A in terms of the 7. Therefore using this property, 
Eq. ( |B.lD can be recast as; 

F{u,uJi\Fig. M) = \ {^y^H j v^j v^ijvj2gco {i)gu>, {ii)gu>-u., (72)^7^71^72 

X5(7l7i,72)r(7|7i,72) (B.3) 
where g (7171,72) represents the integral over the center of mass coordinates of the three BFKL states; 

g (7l7i,72) = jd^nj d^Ri I d^R2iRoiRoif'" {Ri2R\2t' iR2oR*2of'° (B.4) 

where [|] /3oi = -]^-]^{n - ni + n2) + {iy - ivi + 1^2) 

P12 = -]^ + ]^{n + ni+n2) - {iy + ivi+iy2) 

(320 = + ^ {-n - rii + n2) + {iu + iui - iiy2) (B.5) 
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Substituting for g^^ (7) the explicit expression of Eq. ( 2.12| ) and switching to rapidity representation 



(using the above explained inequality y ' + Syn < yi < 2/ for the integration limits of the vertex rapidity 
2/1); 



F{A,6yH\Fig.E) = dyj ^.e^^ |^ e^^^i-.'+^.n^^^ 

= (^Yah r dyi fvjXij) e-(^)^ 

j T^lij ^72A (71) A (72) exp{(u;(7i) +^(72) - ^(7)) (yi - y') -^(71)^^//} 
xE^E^,E^,g (7|7i,72) T (7I71, 72) (B.6) 

where A = y — y' is the rapidity gap between the scattering dipoles in Fig. The leading contribution 
stems from the region where all conformal spins {n, 711,712} = since the BFKL eigenfunction uj{n,u) 
decreases sharply as n increases, and is positive only for ti = at high energies. Therefore ignoring all 
solutions where {n, 711,712} 7^ 0, and using the explicit expression of Eq. ( |2.1[l| ) for the integration measure 
Eq. ( |B.6|) reduces to; 

-i) Ah dyJ duh{u)\{v)e^^-^'' 

4 ^ Jy'+Syn J -00 



Ah 



dv 



287riV3y_^ {l/2 + ivf{l/2-iu)' 



.e'-(-)^E,fiu\A,6yH) 



(B.7) 



where f {v\/S.,5yH] 



1 



y 



dui v\ 



[1/2 + ivif {1/2 + ivif 



dv2V2 



3(1/2 + ii/2) (1/2- ii/2) 



X / dyiey.Y){{uj{ui)+uj{v2) -uj{y))[yi-y') -uj{vi)5yH} 
'y'+SyH 

xEy^Ey^g {h'\h'i, V2) r {u\ui, V2) 



(B.8) 



Eq. ( |B.§| ) represents the contribution of the fan part of the diagram of Fig. |2^. The leading contribution 
to Eq. ( |B.8| ) stems from the region 1^1, 1^2} = {1/2,0,0} which was discussed in the previous sections 



for solving the Pomeron loop integrals. For this region, g{u\{),Qi) can be calculated from Eq. ( B.4 ), by 
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choosing the frame of reference where one of the center of mass coordinates R = 0. With this choice of 
frame after evaluating the remaining integrations over Ri and R2 it turns out that; 



After substituting Eq. (|B.g| ) and the triple Pomeron vertex of Eq. (|A.19f) into Eq. (lB.g|) , the integrals 
over ui and U2 can be solved using the method of steepest descents. Plugging the result back into Eq. (| 
yields the expression; 



X / dyie'^^'^^y-y^^ ^ (B.IO) 



where C is the contour shown in Fig. 13 which encloses the pole at iu = 1/2, and a factor of 2 has been 
included which takes into account the identical contribution from the contour C , which encloses the pole 
at iv = —1/2. One can solve the u integral using the same technique shown in Eq. ( |3.24| ), which generates 
the derivative of the Dirac delta function 5^'^^ (y — yi) /^s, and this is absorbed by the integration over the 
rapidity variable yi to finally yield; 



F{A,6yH\Fig.l 



26? Ah ( e2-(0){A-5s/H/2) 



N^^[u;" (0)]3 dA? (A - 5yHf^ 



2.93 X 10"^^ GeV"^ (a, = 0.12) 

3.95 X 10"^ GeV"2 (q^ = 0.2) (B.ll) 
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